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Abstract 

The first example of the so-called "coupled" integrable hydrodynamic chain is 
presented. Infinitely many commuting flows are derived. Compatibility conditions 
of the first two of them lead to the remarkable Manakov-Santini system. Integrabil- 
ity of this four component three dimensional quasilinear system of the first order as 
well as the coupled hydrodynamic chain is proved by the method of hydrodynamic 
reductions. In comparision with a general case considered by E.V. Ferapontov and 
K.R. Khusnutdinova, in this degenerate case component hydrodynamic reduc- 
tions are parameterized by + M arbitrary functions of a single variable, where M 
is a number of branch points of corresponding Riemann surface. These hydrody- 
namic reductions also are written as symmetric hydrodynamic type systems. New 
classes of particular solutions are found. 
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1 Introduction 



This paper is dedicated to a new type of integrable systems, i.e. vector hydrodynamic 
chains. The simplest example of scalar integrable hydrodynamic chains is the famous 
Benney hydrodynamic chain (see [2]) 



= ^^+1 + kA'^-'Al, A; = 0, 1, 2, 



(1) 



This hydrodynamic chain can be written in the conservative form^ (see, for instance, [31]) 



dtHk = [Hk+i HnM^i-S] , A; = 0, 1, 2, 



(2) 



m=0 



where all conservation law densities Hk are polynomial functions with respect to moments 
A". For instance (see [2]), Hq = Hi = A\ H2 = + {A^, = A^ + 3A°A\ ... 
In this paper, we consider a simplest generalization of the Benney hydrodynamic chain 



to a vector case 
/ \ 



( \ 

y fk+i j 



( \ 

(Ik 



fc-1 



m=0 



1— m 



y fm J 



k = 0,1,2 



, J. , . . . , 



(3) 



where the corresponding vector function contains two species of field variables {qk, fk) 
only. We believe that much more complicated integrable vector hydrodynamic chains 



A^ = 5o(A)i^.+i + fii(A)A^ + ... + Bk+iiA)Al A; = 0, 1, 2, 



(4) 



can be discovered in coming future. In such a general construction, matrix functions 
i?fc(A) can depend on A'^,A^,...,A'' for each index k only, and vector A'^ can depend on 
an arbitrary number of components. 

This vector hydrodynamic chain (3) possesses two obvious reductions. If fk vanish, 
then qk Hk, and (3) reduces to (2); if gk vanish, then (3) reduces to the well-known 
linearly degenerate hydrodynamic chain (see formula (37) in [2G]) 



(/fe)t = (/fc+i). - /o (AO. , fc = 0,l,2. 



(5) 



^as usual, a summation for negative values of upper indices should be ignored in all cases below. 



Thus, (3) is a natural generalization of two remarkable integrable hydro dynamic chains 
on "two species" case. This paper is particularly inspired by the first example in physical 
literature such a Taranov's bi-chain (describing high frequency electron-positron plasma 
waves, see detail in [31]) 

A^^ = Al+^ + kEA^^\ = B^+^ - kEB''~\ k = 0,l, 2, 

where 

dE = (A° - B^)dx + {A^ - B^)dt. 

However, this bi-chain has no a "hydrodynamical" origin, i.e. this bi-chain does not pos- 
sess hydrodynamic reductions (see [11]). In a contrary with this bi-chain, hydrodynamic 
bi-chain (3) possesses infinitely many hydrodynamic reductions. Moreover, hydrodynamic 
bi-chain (3) is connected with the three dimensional quasilinear Manakov-Santini system 
determined by a commutativity condition of two vector fields (see [21]). In this paper, 
these vector fields are utilized for a derivation of hydrodynamic bi-chain (3). Method of 
symmetric hydrodynamic reductions (see [30]) is extended on this hydrodynamic bi-chain. 
A construction of higher symmetries (see [3]) is reformulated for a description of hydro- 
dynamic reductions, higher commuting flows and for an application to the generalized 
(Tsarev's) hodograph method (see [35]). Thus, in this paper, we prove an integrabil- 
ity of hydrodynamic bi-chain (3) by the method of hydrodynamic reductions adopted 
on two species of moments A^ and B^ . Moreover, extracting symmetric hydrodynamic 
reductions (see [30]), we present a list of infinitely many explicit hydrodynamic reduc- 
tions, i.e. one can construct infinitely many particular solutions for the aforementioned 
Manakov-Santini system as well as for hydrodynamic bi-chain (3). First examples of 
three-component hydrodynamic reductions of the Manakov-Santini system were found in 
[6]. In comparison with the previous preliminary investigation (see [(:>]), we been able to 
prove a consistency of all described hydrodynamic reductions with hydrodynamic bi-chain 
(3) as well as with a whole Manakov-Santini hierarchy. 
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This paper is organized in the following way. In Section 2, a four component three 
dimensional hydrodynamic type system (equivalent to the Manakov-Santini system) is 
presented. Due to its dispersionless Lax representation, a first example of integrable 
vector hydrodynamic chains is derived. Simple reductions to well-known hydrodynamic 
chains are found. In Section 3, a generating function of conservation laws and commuting 
flows is described. In Section 4, the Manakov-Santini system is obtained from a new hy- 
drodynamic bi-chain. In Section 5, semi-Hamiltonian reductions of the Manakov-Santini 
system are considered. It is proved that a corresponding dispersion relation is degen- 
erated. In Section 6, two species of Riemann invariants are introduced. Full extended 
Gibbons-Tsarev system is derived. In Section 7, all symmetric hydrodynamic reductions 
are described, and some of them are explicitly found. In Section 8, the generalized hodo- 
graph method is applied for a construction of infinitely many particular solutions for the 
Manakov-Santini hierarchy. In Conclusion, a generalization of the approach presented in 
this paper is discussed. 

2 Lax representation 

The function X{x, t, y; q) determined by the pair of linear equations 

At = (g - a)\^ - u^Xg, \y = {q^ - aq - c)\^ - {qu^ + Ut)\ (6) 

exists, because a compatibility condition {\t)y = {\)t implies to the three dimensional 
hydrodynamic type system 

at = (c + u)x, ct = cttx - acx - + ay, ut = - aux-, h = cu^ + Uy, (7) 

which is equivalent to the well-known Manakov-Santini system (see [21]) 

Wxy = Wtt + WxWxt + {U- Wt)Wxx, Uxy = Utt + ul + WxUxt + {u - Wt)Uxx- (8) 
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Indeed, the first two equations in (7) reduce to the first equation in (8) by the potential 
substitution a = Wx and c = Wt — u, while the second equation in (8) can be obtained 
from the compatibility condition {bt)x = {bx)t of other two equations in (7). 

An inverse transformation X{x,t,y;q) — > q{x,t,y; X) leads to the pair of compatible 
quasilinear equations 

qt = {q- a)qx + u^, qy = {q^ - aq - cjq^ + qUx + Uf (9) 

Let us introduce another compatible pair of auxiliary quasilinear equations [formally, 
replacing qx, qt, qy and u by fx, ft, fy and a, respectively) 

ft = {q- a) fx + ax, fy = {q^ - aq - c)fx + qax + at, (10) 

where f{x,t,y;q) depends on the parameter A implicitly via a dependence q{x,t,y; X). 
The compatibility conditions {qt)y = {qy)t and {ft)y = {fy)t imply to (7) in both cases. 

The reduction a = 0, c = —u reduces (7) to the famous Khokhlov-Zabolotskaya equa- 
tion (also well known as a dispersionless limit of the Kadomtsev-Petviashvili equation) 

ut = hx, ht = Uy-uUx (11) 

determined by compatibility condition (9) under the reduction / = 0; the reduction 
M = 0, 6 = reduces (7) to the remarkable equation (recognized as an anti-self- duality 
reduction of the famous Yang-Mills equation, see [8]) 

at = Cx, Ct = cttx - acx + ay, (12) 

determined by compatibility condition (10) under the reduction q = 0. 

Let us substitute an expansion (A oo,q oo, according to a general theory of 
a relationship between dKP equation (11) and Benney hydrodynamic chain (2); see, for 
instance, [31]) 



\ ^0 ^1 ^2 
q = X-j-^-^-... (13) 
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in the first equation of (9). Tlien tlie first part of hydrodynamic bi-chain (3) is derived, 
wliere u = qo and b = qi. A substitution of an expansion (A oo, / oo, according 
to a general theory of a relationship between quasilinear system (12) and hydrodynamic 
chain (5); see formula (39) in [26]) 



in the first equation of (10) implies to the second part of hydrodynamic bi-chain (3), 
where a = fo and c = /i — go — /o/2. Thus, whole hydrodynamic bi-chain (3) cannot be 
derived just from the first equation of (9) as well as from the first linear equation of (6) 
by an inverse transformation (A — oo, q — oo, see (13)) 

^° 

X = q+ — + — + —..., (15) 
q q^ q-^ 

which is useful in a derivation of Benney hydrodynamic chain (1) under a corresponding 
reduction / = to the dKP case (see, for instance, [31]). In such a case, another formal 
expansion (A — oo, M — >^ oo, the so-called Orlov function; see [3]) 

oo oo 

M = Y^ ne~^r~^ + w„A-" 

n=l n=X 

for the first equation^ of (6) can be utilized in a derivation of the second part (see [6]) 
of hydrodynamic bi-chain (3), where p„ = {vn)x are nothing else but conservation law 
densities (see below). Here t° = x,t^ = t, = y. All other higher "times" are associated 
with higher commuting hydrodynamic bi-chains (see below). Instead this Orlov function 
M, in this paper, we use an alternative approach (see (9) and (10)) to derive a generating 
function of conservation laws {p = M^, where A is a free parameter; see detail in [6]) as 
well as a generating function of commuting flows (see below). 

A substitution (15) to (14) implies to the expansion (/ — ^ oo, g — > oo) 

f = 2~~~- (16) 

q g-* 



^It means that the function M(A,t,v) satisfies (6) as well as the function A(g, A). 















( \ 
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Let us rewrite (10) splitting dynamics with respect to x,t,y and q. It means the 
dynamics (9) is excluded from (10), and q becomes an extra independent variable as well 
as x,t,y. Corresponding equations reduce to the form (cf. (6)) 

ft = {q- a) fx - u^fg + a^, fy = {q^ - aq - c)fx - {qu^^ + Ut)fq + qa^ + at. (17) 

A substitution (15) to the first equation in (6) and (16) to the first equation in (17) 
leads to a generalization of Benney hydrodynamic chain (1) to the vector case (cf. (3)) 



y4^, k — 0,1, 2, 



where a = and u = A^. 

This hydrodynamic bi-chain has three obvious reductions. If A'' = 0, then (18) reduces 
to (5); if B'' = 0, then (18) reduces to (1); if B'' = A'', then (18) reduces to the "deformed" 
Benney hydrodynamic chain 

— ^x ^ ^x ^ ^xi ft- — U, i, Z, 

derived in [13], whose integrability was investigated in [32]. This hydrodynamic chain is 
connected (see [29]) with the so-called "interpolating system" (see [7]). 

Moreover, Manakov-Santini system (8) possesses a natural reduction u = Wt (see [6]) 
to the continuum limit of the discrete KP hierarchy (see, for instance, [20]) 
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where w = Vlf . In such (7) reduces (i.e. c = 0) to 

o-t = dy = bx, bt = Uy, Ut = bx — aUx- 



(19) 



A comparison of the first conservation law at = Ux (i.e. B^ 
(18)) with the first conservation law of (18) 



y4°, see the notation in 



B't = - liB^) 



implies to the first constraint 



It means, that hydrodynamic bi-chain (18) reduces to the hydrodynamic chain 

= - + kB''-^ (^B^ - ^{By^ , (20) 

which is a particular case of a more general class of integrable hydrodynamic chains (a, /9, 7 
are arbitrary constants) 

= 5^+1 - + {ak + f3)B''Bl + ^kB^"-^ (^B^ + ^~^~^ (E°)^^ 

investigated in [18] and then in [32]. A computation of all other constraints A^{B^, B^, B^ 
is a more complicated task. Nevertheless, the second constraint can be found due to 
comparison of the fourth equation Ut = — aux (see (19)) with the first equation 

= A], — B^A^ from the A-part of (18). It means that b = A^ and a substitution 
of the first constraint A^ = B^ — (5^)^/2 leads to the conservation law B^ = A].. Taking 
into account the second equation from (20), the second constraint A^ = B^ — (5*^)^/3 can 
be easily obtained. 

In aforementioned case (19), the first equation in (9) reduces to 

Pt = dxieP + B'), 

where p = ln(g— 5°). This is a generating function of conservation laws for the remarkable 
hydrodynamic chain (see, for instance, [15] and [20]) 

= + A:C"=C°, A; = 0,1, 2,..., 

where C° = B^ and all other moments are connected with 5" by inverse triangular 
point transformations (see detail in [32]). Since this hydrodynamic chain belongs to the 
Boyer-Finley (a continuum limit of 2D Toda Lattice) hierarchy (see [5]), a corresponding 
algebraic mapping is given by the expansion {X ^ oo,q ^ oo,p ^ 00) 



which must coincide with two other expansions (A — >■ oo, q — >■ oo; see (15)) 

A = g + — + — + — ... = (g - 5°) exp — + — + — ... . (21) 
q q-i \ q q'^ q-^ J 

Thus, we see that hydro dynamic chain (20) is embedded in hydrodynamic bi-chain (18) 
by very comphcated reductions A^{B^ ,B^, ...^B^'^^), while (8) reduces to (19) by the 
simple constraint c = (see (7)). Moreover, in this case, both linear systems (9) and (17) 
coincide due to aforementioned link (see (16) and (21)) 

A = (g - 5°)e-^. 

The general problem of all possible compatible reductions A^{B^ ^ B^i 5^+") should 
be considered in a separate publication (for each fixed non-negative integer n). 



3 Generating functions 

Let us rewrite the first equation of (10) in the form 

st = {q- a)s^ + sa^ 

where (see (14)) 



, = exp/ = l-|-il-£|-... (22) 



Then hydrodynamic bi-chain (3) reduces to the form 

k-l k-1 

{Sk)t = iSk+l)x — So{Sk)x — ''^^qk~l-m{Sm)x + Sk{So)x, {Qk)t = iqk+l)x — So{qk)x — ''^^qk~l~miqm)x, 

m=0 m=0 

where Sq = /o (all other higher unknown functions Sk are polynomial functions of fo,qo, fi, 
qi, fk,qk for each index k). If qk vanish, this hydrodynamic bi-chain reduces precisely 
to the form derived in [1]. 

Let us introduce the so-called vertex operator {( oo) 

dric) = -^dto - ^dti - ^dt2 - ... (23) 
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Theorem: Higher commuting hydrodynamic bi- chains are determined by the pair of 
vertex equations 

driO^W = '^^^^^'/(A) ~ 'j^f ^'^^^ ^ d^iO^W = s{C)d. ln(g(A) - g(C)) (24) 
by a substitution (13), (22), 



c c 



and (23). 

Proof: A direct substitution (13), (22) and (23) in (24) leads to a family of in- 
finitely many hydrodynamic bi-chains. Their commutativity follows from the compatibil- 
ity conditions dr{Q{dtS{X)) = dt{dr((;)S{X)), dr{Q{dyS{X)) = dy{dr{(;)S{X)), dr{Q{dtq{X)) = 

dt{dT-[()q{X)),dr[(){dyq{X)) = dy{dT-[()q{X)) leading to the extra vertex dynamics (see (7)) 
dr{c)a = d^s{C), dr{QC = (g(C) - a)d^s{C) - s(C)9^(g(C) - a), 

(25) 

fq\n \ 

dr(c)U = s{C)d^q{C), dr{c)b = s{C)d^ ( + n 1 , 

which are compatible (i.e. dr{(;){dta) = dt{dr{(;)a), dr{Q{dtb) = dt{dr{Qb),dr{Q{dtc) = 
dt{dr(()c),dr((){dtu) = dt{dr(()u)) with (7). Thus, the theorem is proved. 

Remark: This generating function of commuting flows (25) can be obtained by a 
substitution (13), (22) in (24) only. 

Taking into account the flrst equation in (25), expansions (22) and (23), the sole 
function ^(t°, t\ ...) can be introduced such that 

dS 

s(C) = drioS +1, Sk = 

Then the flrst equation in (24) implies to^ 

s{X)d.J{X) - s\X)d^s{X) 



driX)S{X) 



^here and below, the sign "prime" means a derivative with respect to a corresponding parameter. 
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in the limit ( ^ X. Thus, all functions qk can be expressed via first and second derivatives 
of the sole function 5* with respect to corresponding "time" variables t". A substitution 
of these dependencies to the second equation in (24) leads to an infinite set of three 
dimensional quasilinear equations of the third order. From the other hand, three copies 
of the first equation in (24) imply to a sole vertex equation 

^ [dr(X)S + l]d:^dr{n)S - [dr{r,)S + l]d^dr(x)S _ Q 

dT{X)dT{ri)S 

An expansion of this generating function with respect to three parameters A, (, rj at infinity 
leads to an infinite set of three dimensional quasilinear equations of the second order. 
Remark: The second equation in (24) implies to 

.3^x _ drix)qiX) 
'^^^ a.lng'(A)- 

in the limit ( \. Thus, all functions Sk can be expressed via functions go, qi, qk and 
their first derivatives with respect to corresponding "time" variables t". A substitution of 
these dependencies to the first equation in (24) leads to an infinite set of three dimensional 
quasilinear equations of the second order. 

\n(\\-n(n]n ^^Wg(^) _ ^^(C)g(C) ^ '9r(A)g(A) gr(A)g(A) gr(c)g(C) 
giUJ^.(C)^^l^^,(^) 9.1ng'(C)'^^a.lng'(A) dMWdJnq'iCY 

Generating functions (24) have an obvious two reductions. The choice s(A) = 1 means 
reduction to the generating function associated with the dKP hierarchy (11) and with the 
Benney hydrodynamic chain (2); the choice g(A) = A means reduction to the generating 
function associated with hydrodynamic chain (5), quasilinear system (12) and the so-called 
"universal" hierarchy (see detail in [1] and [26]). 

Theorem: Hydrodynamic bi-chain (3) possesses a generating function of conservation 
laws 

dtp = d:,[{q - a)p], (26) 
12 



where the generating function of conservation law densities is given by 

P(A) ^ fl. (27) 

Proof: Indeed, taking into account the first equation in (9) and the first equation in 
(10), a direct substitution (27) in (26) imphes to an identity. 

Remark: A substitution (13) and (22) in (27) leads to an expansion 



where all coefficients pk can be expressed polynomially via field variables qo, fo,qi, fi, ...,qk, fk- 
Corresponding conservation laws can be written in the form 

Pk,t = dx I Pk+i - PoPk -Qk- ^Pfc-n.-ign I , /c = 0, 1, 2, ... 

\ n=0 / 

Of course, all conservation laws can be written via moments A'', B^. For instance, first 
three of them are 

d,B' = dx{B'-^-^^, d,[B'+^^ + A'^=dx{B--^^ + A^ 
dtlB^ + B'B' + + 2A' + 2A'B^ 

djB' + ^ - _ W ^ + 2A^B' - {B^A^ + {A^ 
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Corollary: A generating function of conservation laws and commuting flows is given 
by an auxiliary vertex equation 

Indeed, a substitution (24) and (27) in the above equation implies to an identity. More- 
over, this generating function can be written via the sole function S 

dx[\ns{X)]' _ [s{X)dxs'iX)-s'iX)dxs{XMC)driosiX) 
drix) In s{X) ^ [s{Odxs{X) - s{X)dxs{CMX)drix)s{X) ' 

or via infinitely many field variables qk 

dxq'jX) g'(C)5xg'(A)9.(c)g(C) 



drWQW [q{X) - q{0]dxq'{C)drwq{X) ' 
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4 Derivation of the Manakov— Santini system 



A substitution of expansions (15) and (16) in the second equations of (6) and (17) implies 
to the first commuting fiow of (18) 



\ / y 



I *^^2 \ ( /lfc+1 \ / /^nxo s / Ak\ 



B^-i^-A» 



A'- 













\ 






+ ik + l) 








{Al-B'Al), fc = 0,l,2,..., 


(29) 










r>k-l 


) 







where c = B^ - {B^f/2 - A°. 

Let us take first two equations from (18) and the first equation from (29) 





A^ 






( A^ 




-B' 


/ 


A^ 


\ 








--{ 












v 




) 






I 


X 


v 


B' 


/ 



B' 



A' 



B^-^^-A" 



An ehmination A^ and B^ leads to three dimensional four component quasilinear system 
of the first order 





A^ 






/ 


A' 






/ 


A^ 


\ 






















v 


B' 


) 


t 


v 


B' 


J 


X 


\ 


B' 


/ 



^A'\ 



B^-i^-A' 



A substitution A° = u, B° = a, A^ = b and B^ = c + a^/2 + u implies to (7). 

Thus, we have proved that three dimensional four component quasilinear system of 
the first order (7) is equivalent to commuting hydrodynamic bi-chains (18) and (29). 
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It means that any solution (obtained by the method of hydrodynamic reductions, see 
below) of (7) can be utilized for a construction of a corresponding solution of two linear 
systems (6) and (10); then such a solution can be used for construction of a solution 
of commuting hydrodynamic bi-chains (18) and (29). And vice versa, any solution of 
commuting hydrodynamic bi-chains (18) and (29) implies immediately to a solution of 
(7). 

5 Semi-Hamiltonian hydrodynamic reductions 

Following approach established in [11], all field variables a,b,c,u in (7) are considered 
as functions of Riemann invariants A^(x, t,?/), where hydrodynamic reductions are 
commuting semi-Hamiltonian hydrodynamic type systems 

Xi = v\X)Xi, y=w\X)Xl, i = l,2,...,N (30) 



for any A^. Characteristic velocities must satisfy the systems (see [3-^)]) 

dkv"- dkW 



du4^ = d,4^, ^7^J^k; (31) 



i ^ k, (32) 



where dk = d/d\^ . A substitution (30) in (7) implies to the so-called "dispersion relation" 

w' = +av' -c (33) 

and two differential relations 

di{c + u) = v'^ditt, dih = (f* + a)diU. (34) 
A substitution (33) in (32) together with the first equation from (34) implies to 

V = -^, ^^k, (35) 
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where we introduce = f * + a instead f*. Compatibility conditions dk{di{c + u)) 
di{dk{c + u)) and dk{dib) = di{dkb) for each pair of distinct indices lead to 

„ diU ■ dkU 



dikU 



k\2 ■ 



(36) 



diktt 



diU ■ dktt + dkU ■ dia 



(37) 



(g* — q'^y 

System (35)-(37) is an analogue of the Gibbons-Tsarev system for the dKP equation 
Equations (37) is a linear system whose variable coefficients are determined by 



[see 

solutions of sub-system (35), (36) which is precisely the aforementioned Gibbons-Tsarev 
system. Original Gibbons-Tsarev system (35), (36) has a general solution parameterized 
by arbitrary functions of a single variable. Linear system (37) also has a general solution 
parameterized by N arbitrary functions of a single variable. Thus, whole system (35)-(37) 
has a general solution parameterized by 2N arbitrary functions of a single variable. 

Thus, this is a first example (in a literature), where number of arbitrary functions 
twice bigger than in an ordinary theory (see [11]). Let us give an explanation of this 
deviation here. Three dimensional four component hydrodynamic type system (7) should 
be written in the form 
/ 



\u J 



1 

c 

c —1 —a 

\0 I -c 



1 \ fa\ 



I 





1 

10 

y y 



{31 



\u J 



A dispersion relation (see [11]) for any three dimensional hydrodynamic type system 



Ut = Aux + Bu^ 



is determined by 



det[i + ^5 -/iE] = 0, 
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where ^ and /i are constants, and E is an identity matrix. In our case, 



det 



-/i 





1 1 







c + ^ 




-1 


— /i — a 





1 


— /U — a 







(39) 



implies to the reducible quartic 



{fi' + afi-c-O' = 0. 



Thus, instead a full algebraic equation of the fourth order (with distinct roots), hydrody- 
namic type system (7) possesses the degenerated dispersion relation (conic) (cf. (33)) 

^ = 11^ + afi — c. 

An existence of hydrodynamic reductions (30) (integrability conditions (31) and (32) are 
not necessary!) reduces (38) to the form (cf. (39)) 

— U \J ± ± 

b 







V 








1 


1 










c + 


c + 


-1 


—v^ — a 








1 





— f * — a 



J 



Thus, indeed (see (33)) ^ has a meaning of the characteristic velocity w, and fi has a 
meaning of the characteristic velocity v. Such an existence of this double conic leads to a 
double freedom of solutions for generalized Gibbons-Tsarev system (35)-(37). 

An existence of such a double conic leads to the following phenomenon. The number 
of Riemann invariants A'^ can be split on two parts. First M Riemann invariants are 
branch points of the algebraic curve A = A(A; q) (see below); all other N — M Riemann 
invariants are just mark points r" on this algebraic curve A = A(A; q). M can run all values 
from 1 up to A^. Another origin of this phenomenon is connected with an observation that 
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whole hydrodynamic bi-chain (18) contains two parts, where each of them is determined 
by a different hnear system (see (6) and (10)). Moreover, first M Riemann invariants are 
connected with the A-part of (18), while other N — M Riemann invariants are connected 
with the 5-part of (18). 

Since integrable three dimensional four component hydrodynamic type system (7) 
determined by linear system (6) is associated with commuting hydrodynamic bi-chains 
(18) and (29), to avoid a lot of unnecessary computations without loss of generality and 
for simplicity we shall consider just first hydrodynamic bi-chain (18) below. 

6 Two species of Riemann invariants 

Following the approach presented in [9], let us suppose that all moments are functions 
of M Riemann invariants A". In such a case, the first part of (18) reduces to an infinite 
series of recursive relations 

d,A^+^ = q'd.A^ - kA^-^d.A'', = 0, 1, 2, 

whose compatibility conditions di{djA^) = dj{diA'^) are nothing else but the original 
Gibbons-Tsarev system (35), (36). Moreover, these Riemann invariants A* are nothing 
else but branch points A* of the Riemann surface A = A(A; q), i.e. A* = A(A; g'), where 
g*(A) are solutions of an algebraic system (see [9] and [22]) 

9,A|,=,. = 0, i = l,2,...,N. 

Indeed, linear system (6) reduces to the so-called Lowner equation (see [9]) 

3-11 

d^A = ^-d,A, (40) 
g — 

whose compatibility conditions di{dkA) = dk{diA) imply to original Gibbons-Tsarev sys- 
tem (35), (36). Thus, hydrodynamic bi-chain (18) reduces to a composition of the inte- 
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grable hydrodynamic chain 

= B^+^ - B^B^ + kB''~\^, A; = 0, 1, 2, ... (41) 

and the semi-Hamiltonian hydrodynamic type system (see the first equation in (30)) 

Aj = (g^(A) - i?0)A^, z = l,2,...,M, (42) 

where functions u and q^{X) satisfy original Gibbons-Tsarev system (35), (36). Never- 
theless, linear system (37) cannot be derived from (6), because 5° is not a function of 
Riemann invariants A'. However, if we consider hydrodynamic reductions of hydrody- 
namic chain (41), i.e. if we suppose that is a function of Riemann invariants A*, then 
all higher moments i?^' must be functions of Riemann invariants too. Then compatibility 
conditions di{djB'') = dj{diB'^) lead to (37), where 

diB^+^ = q'diB^ - kB^^^diU. 

A most interesting consequence of aforementioned degeneracy described in the previous 
Section is an existence of extra Riemann invariants which are just mark points on already 
determined Riemann surface A = A(A; q). Indeed, let us introduce N — M mark points 
= A(A; q^) such that each inverse function q^{X, r^) is a solution of the Lowner equation 
(see (40), where we consider an inverse function g(A; A)) 

3'ZL 

d,q=-^, z = l,2,...,M. (43) 

q'-q 

A general solution of this overdetermined system (i.e. the compatibility conditions di{dkq^) = 
dk{diq^) must be fulfilled due to original Gibbons-Tsarev system (35), (36)) 

d,q^ = -^, z = l,2,...,M, j = l,2,...,iV-M (44) 
— q^ 

depends on M arbitrary functions of a single variable (because this is a system of the first 
order in partial derivatives with respect to M independent variables A*) and other N — M 
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arbitrary functions of a single variable (because each Riemann invariant r'' as a mark 
point on the Riemann surface A = A(A; q) is determined up to an arbitrary transformation 
— > Rk{r'^). Let us remind that a solution A of linear system (6) is determined up to 
an arbitrary transformation A A(A)). Then let us introduce the functions bj{X,r^) 
determined by their first derivatives (the compatibility conditions di{dkbj) = dk{dibj) are 
fulfilled due to original Gibbons-Tsarev system (35), (36)) 

such that 

B°{\r)=Y^ expcj{X,r^)dr^ (46) 

satisfies overdetermined system (37). Whole overdetermined system (37), (44), (45) can 
be rewritten as a quasilinear system of the first order on unknown functions , diB^ and 
djB^ = expcj (here and below in this Section dj = d/dr^). Thus, its general solution is 
parameterized by arbitrary functions of a single variable. It means, a general solution 
of complete overdetermined system (35)-(37), (44), (45) is parameterized by + M 
arbitrary functions of a single variable. 

Let us introduce component hydrodynamic type system (see (42)) 

K = {q\\)-B')\l, ri = {q^{\r^)-Byi. (47) 

An integrability condition (i.e. semi-Hamiltonian property (31)) is fulfilled automatically 
due to (35)-(37), (44), (45). 

Theorem: Hydrodynamic bi-chain (18) possesses infinitely many N component hydro- 
dynamic reductions (47) determined by the extended Gibbons-Tsarev system (35)-(37), 
(44), (45). Then all higher moments B'^ can be reconstructed iteratively in quadratures 

M N~M 

dB^^^ = "^[q'diB^ - kB^~^diu]dX + ^ q^djB'^dr^ A; = 0, 1, 2, ... (48) 
i=i j=i 
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Proof: Since all moments depend simultaneously on both species of Riemann 
invariants A* and H, then the 5-part of (18) reduces to the second equation in (47), 
where (see (48)) 

d,B''+^ = q'diB^ - kB^~^diU, djB''+^ = q^djB''. (49) 

The compatibility conditions (9j(g*9ii?'') = di{q^djB^),di{q^djB^) = dj{q^diB^—kB^~^diU), 
dj{q^diB^ — kB'^~^diu) = di{q^djB^ — kB^^^dju) lead to the system [i ^ j in the first and 
third equations below) 

didjB'' = 0, {q' -q^)d,djB^ = diq^ ■djB'' + kdiU-djB^-\ 

(50) 

(g* - q^fdidjB'' = diu ■ djB^ + dju ■ diB'' + k{q' - q^){diU ■ djB^-^ - djU ■ diB''"^). 

The second equation in (49) can be written in the common form 

djB'' = {q^fdjB^, A; = 0, 1, 2, ... (51) 

A substitution (51) in the first two equations of (50) leads to (44), (45) and (46), where 
each function q^ depends on all Riemann invariants A* and just one Riemann invariant rK 
The third equation in (50) coincides with (37) for A; = 0. The compatibility conditions 
di{djiB^) = dj{diiB^) can be verified by the induction principle (taking into account the 
first equation in (49). Precisely, such a computation was made for A-part of hydrodynamic 
bi-chain (18), i.e. for Benney hydrodynamic chain (1) in [9]. 

Correspondingly, linear problem (10) determining the -B-parts of hydrodynamic bi- 
chains (18) and (29) reduces to (cf. (40)) 

a,F = 4^, (52) 
q — q Q-' ~ Q 

where F{X, r) is used instead of f{x, t, y, q) for hydrodynamic reductions (30). The depen- 
dence g-'(A, r-') and (46) follow from the compatibility condition di{djF) = dj{diF). The 
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compatibility condition di{djF) = dj{diF) leads to (44) and (45). The compatibility con- 
dition di{djF) = dj{diF) satisfies automatically due to original Gibbons-Tsarev system 
(35), (36). Similar computations can be repeated for generation function of conservation 
laws (26). Then the generating function of conservation law densities p can be found in 
quadratures 

M / r^ o \ N-M X 

j-l\{q-q'Y q-q'J j^i q - q' 

It is well known that generating function of conservation laws (26) is associated with 

the so-called "linearly degenerate" hydrodynamic chains and their linearly degenerate 

hydrodynamic reductions and (see [26] and Section 9 in [30]). First such an example 

was the Whitham equations of the Korteweg de Vries equation (see [16]). More general 

theory was suggested in [17]. Theory of linearly degenerate hydrodynamic type systems 

is presented in [10], [12], [25]. Following [16] and [17], let us introduce the so-called 

"quasi-momentum" differential dP = pdX and the so-called "quasi-energy" differential 

dQ = {q — a)pdX. Then characteristic velocities are given by 



where the differentials of "quasi-energy" and "quasi-momentum" possess similar singular- 
ities on the Riemann surface A = A(A; q) at vicinities of mark points A = r*. 

Generating function of conservation laws and commuting flows (28) can be written 
in similar form. In such a case, the generating function of "quasi-energy" differentials is 
given by 

dQ{X,0 = W{XX)p{X)dX ^ -^^^^^p{X)dX. 

qW - q{0 

Corresponding generating function of commuting hydrodynamic type systems (see (47)) 
is given by (replacing g(A) in the expression for W{X, Q by and , respectively) 



, _ expF(C) , _ expF(C) 
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This hydrodynamic type system possesses generating function (28). In such a case (cf. 
Section 9 in [30]), 

M r. jjr N-M ~ 

i=l j=l 

must coincide with (53). Indeed, a direct substitution of expressions W{X, Q, W''{X; Q and 
W^{X, r; () implies to (53). Moreover, hydrodynamic type systems (47) and (54) commute 
to each other (see Tsarev's condition (32)). It means, that any hydrodynamic reduction 
(47) satisfying extended Gibbons-Tsarev system (35)-(37), (44), (45) is compatible with 
the whole Manakov-Santini hierarchy described in Section 3. 

Remark: Semi-Hamiltonian property (31) and commutativity condition (32) for hy- 
drodynamic type systems which are hydrodynamic reductions of linearly degenerate hy- 
drodynamic chains can be simplified (see Section 9 in [30]) to the form, respectively 



where dP = p{X)dX is a differential of "quasi-momentum" , dQ = V{X)p{X)dX is a differen- 
tial of "quasi-energy" and dQ = W{X, Qp{X)dX is a generating function of "quasi-energy" 
differentials. In the case of two species of Riemann invariants (see, for instance, (47) and 
(54)), the above formulas must be written in the form 

^ dkV ^ d,V ~ dkV a d,V .,, ~ dkV . d,V 
di— — — = dk—. — —, di:~- — = dk:--. — 77, ^ r k; di— — — = dk:~, 



Vk-v '"yi-V 'vk_v yi-V' V^-V v^-v' 

dkV _ dkW dkV _ dkW 

yk-v W^-W V^-V W'^-W' 

7 Symmetric hydrodynamic reductions 

Benney hydrodynamic chain (1) possesses the special hydrodynamic reduction (see [32]) 



A" 



X 
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where all moments are determined by the choice (e, are constants) 

^ M M 

^' = ;fc^E^^(«')'^'' ^ = 0,1,2,...; J]q = 0. 

1=1 1=1 

Let us extend this decomposition on both species of moments (e^ and 7^ are constants) 

M , N-M 



^' = ^E^^(«')'^'' = ^ E ^■(^■)'^'' fc = 0,l,2,... (55) 

i=l j=l 

Then bi-chain (18) reduces to the hydrodynamic type system 

al = {a^ - B^)al + bl = {V - B'H + A", (56) 

where 

M N-M 

In a contrary with semi-Hamiltonian hydrodynamic reductions (see (41) and (42)) of 
the A-part of bi-chain (18), the moment decomposition (see (55)) 

N-M 

i=i 

for the -B-part of (18) leads to the hydrodynamic chain 

^ ^k+i _ ^o^k ^ kA'^-iAl = 0, 1, 2, ... 

equipped by the non-diagonalizable (i.e. non semi-Hamiltonian!) hydrodynamic type 
system (cf. (56)) 

bl = {V -B')b>, + Al J = 1,2,..., N-M. (57) 

Indeed, Riemann invariants r'^(b) cannot exist, because the last term cannot be elim- 
inated by any point transformations. Moreover, in a general case, the A-part of (18) 
reduces to the M component hydrodynamic type system 

ai={a'-B'>)ai + Al z = l,2,...,M, (58) 
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where all moments A'' depend on field variables a" only, and the function u = A'^ satisfies 
the original Gibbons-Tsarev system (see [30]) 

(a* — ci^)uik + UiSuk — UkSui = 0, i ^ k, (59) 

where 6 = Hd/da"^ is a shift operator, Ui = du/da^, Uik = d'^u/da^da'^ , all higher moments 
A'' can be reconstructed iteratively (here and below diA'' = dA^ /da^) in quadratures 

M 

dA^+^ = J^l'^'diA'' + {6A'' - kA''-^)diAya\ A; = 0, 1, 2, ... 
Linear system (6) reduces to the Lowner equation (cf. (40); more detail in [30]) 

aA + ^fl + f ^)"8,A = 0, (60) 

\ m=l / 

written via field variables a^, where the equation of Riemann surface A = A(a;g). The 
compatibility conditions di{djA) = dj{dih) as well as the compatibility conditions di[a^ djA''+ 
{5 A'' - kA^-^)djA^] = djla'diA'' + {6 A'' - kA''-^)diA^] lead to (59). 

In a contrary, the i?-part of (18) reduces to (57) if all moments depend on field 
variables 6" only. It means, that in a general case, moments B'^ must depend on two 
species of field variables a" and b"^ simultaneously (as well as on two species of Riemann 
invariants A* and r^). Just in such a case, hydrodynamic type system (57), (58) can be 
semi-Hamiltonian. Indeed, the i?-part of (18) 

B^ = B^+^ - B^B^ + kB^'-^Al, A; = 0, 1, 2, ... 

reduces to (57) and (58), if (here and below di = d/db\S = T.d/db"^) 

S.fifc+i = a'd.B'' + {6B'' + 5B^ - kB^~^)diA\ djB^+^ = Vd^B^. (61) 

The last equation can be written in the common form (cf. (51)) 

djB^ = {VfdjB^, A; = 0,1, 2,... 
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The compatibility conditions dj[a^diB^ + {6B^ + 6B^ - kB^-^)diA^] = dillPdjB''], 
dj[a^diB^ + {5B^ + ~5B^ - kB^-^)diA^] = di[a^djB^ + {5B^ + 5B^ - kB''-^)djA^] and 
di{VdjB^) = dj{UdiB'^) lead to the system (cf. (50); i 7^ j in the first and third equations 
below) 

d.d^B^ = 0, (a* - b>)didjB'^ + ■ 4(55° + ~5B^) = 0, (62) 
[a' - a^)didjB^ + diA^ ■ dj{6B'' + ~dB'' - kB''~^) - djA'^ ■ di{5B^ + ~5B^ - kB''"^) = 0. 
The second equation in (62) reduces to 



TV -A/ 

(64) 



dfx, + (l + y ~df\ = 0, (63) 

\ m=l / 

which can be also obtained from (60) formally replacing A on — M functions \j and 
correspondingly q on N — M field variables V . In such solution of the first 

equation in (62) is given by (cf. (46)) 

A dependence of 5° with respect to field variables is given by the last equation in (62) 
{a' - a')d,djB'' + ■ 9^(55° + 55°) - 9,A° ■ 9,(55° + 55°) = 0, i ^ j, (65) 
while all other higher expressions 5^(a, b) can be found iteratively (see (61)) 

M N-M 

dB^+^ = ^[a^9,5'= + (55^= + 55'^ - kB''-^)diAya' + ^ VdjB'^dV, k = 0,l, 2, ... 
i=i j=i 

The Lowner equation (60) under the inverse transformation A (a; q) — > g(a; A) (cf. (40) 
and (43)) reduces to the form 

a — q \ , a"^ — q 

\ m=l / 

and the generating function of conservation law densities p (see (26)) can be found in 
quadratures (cf. (53)) 

dlnp = ^ dm{q — a) ■ d\n{q — a"^) — ^ 9„a ■ d\n{q — 6"), 

m=l n=l 
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while (cf. (52)) the function F(a, b) is given by another quadrature (see (22) and (27)) 

M N-M 

dF = Y^ dma ■ dln{q - a"^) + ^ 9„a ■ dln{q - 6"). 

m=l n=l 

Remark: System (65) automatically satisfies (for any solution A^{a) of system (59)) 

if 

M N-M 

B' = J2lma"'+ Y.Pnh\ (67) 

■m=l n=l 

where 7^ and /3„ are arbitrary constants. In such 

M N-M M N-M 

F=J2^m Hq -an+Y. ln(g - 6"), P = 3^ n - «™)"''" n - ^" 

m=l n=l ^771=1 n=l 

All moments A'^' can be expressed via field variables a* as well as Riemann invariants 
A-'. Thus, each field variable a* is a function of all Riemann invariants A", and vice versa 
each Riemann invariant A* is a function of all Riemann invariants aJ^. Since all moments 

are functions of both species of field variables a* and V as well as both species of 
Riemann invariants A* and r-', finally, we need to find a dependence of field variables 
via Riemann invariants A* and rK Taking into account (47), (57) and (58) imply to 

(g^' - 6") ■ 9^-6" = 0, j,n = 1,2,...,A^-M (68) 

and {di = d/dX^,dj = d/dr^ in (68) and below in this Section) 

a.a'^ = -^, d,b- = ^^, ^,k = l,2,...,M, 

which can be obtained directly from (43) replacing q by and 6", respectively. Relations 
(68) imply to the simple dependence = q^{X,r^). Since each Riemann invariant r* is 
a mark point on the Riemann surface A = A(A; q), i.e. = A(A; q^), we conclude that 
equations (58) are connected with (42), while equations (57) reduce to the diagonal form 

rl = {V - a)ri, 

where Riemann invariants = A(A; b^),j = 1, 2, — M. Thus (see (63)), Aj(a; IP) = 
Rj{r^), where Rj{r^) are arbitrary functions, because Riemann invariants are determined 
up to an arbitrary transformation — s> Rj{r^). 
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In the next Section, explicit hydrodynamic reductions are considered. A relationship 
between field variables a'',b'^ and Riemann invariants X^,r^ is investigated in detail. 



8 Explicit hydrodynamic reductions 

A general case contains an arbitrary number of field variables of two species and 6" 
in component hydrodynamic type system (56). The most general reduction of X{p) 
known recently (see [4], [17], [14]) is a combination of rational and logarithmic functions 
with respect to p, i.e. 



K 



k=l 



(1) 



M 

,(3) 



5fc In (^p 



m=l 



where and (5„ are constants, a^""* are functions. The approach allowing to extract more 
complicated reductions is presented in [33]. In this Section a simplest case is considered. 

A^ parametric solution of original Gibbons-Tsarev system (59) given by (see [32]; here 
all constants ej are independent, no such a constraint Se^ = in general) 

M 



J2 (69) 



u = 

m=l 

possesses to reconstruct an equation of the Riemann surface A = A(a; q). A substitution 

(69) in (60) leads to the so-called waterbag reduction (see [5], [9] and [32]) 

M 

A = g- J]e„,ln(g-a'"). (70) 

m=l 

The Riemann invariants A* (a) are branch points of the Riemann surface determined by 

(70) , i.e. the condition (Ag = 0) 

M 



g" — a' 

m=l 
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leads to M expressions ^^(a) as well as to M inverse expressions a*(q). Corresponding 
Riemann invariants A' (a) are given by (70) 

M 

V = g*-5^e^ln(g*-a'"). (71) 

m=l 

Other N—M Riemann invariants r-' (a, b) (see (63)) are mark points on the aforementioned 
Riemann surface 



M 



\\q=w = V - ^e„ln(6- 



m=l 



Since \j are arbitrary functions of corresponding Riemann invariants (i.e. \j = Qj{r^)), 
the function a (see (65)) can be found from the system 

did^a H — — — :d^5a — — :di5a = e^ej > -—. ; — ; — r, i ^ ?, (72) 

m=l ^ ^ 

where (see (64)) 

dja = A, = Q, -Y^em HV - a™) ] . (73) 

\ m=l / 

To avoid a complexity of further computations (an integration of arbitrary functions 
Qj{r^) in (73)), let us restrict our consideration on a first nontrivial case Qj{r^) = l3jr\ 
where Pj are arbitrary constants. Then (73) leads to 

^ N-M M N-M 

a = 5(a) + Pnibl' + E E '"^^nia"- ' ^""^M^^ ' ' ^l' (^4) 

n=l 771=1 rt=l 

where the function -B(a) satisfies (see (72)) 

This linear system possesses a general solution parameterized by N—M arbitrary functions 
of a single variable. Its first nontrivial solution is given by 

M M 

2 



m=l m=l 



where 5 and 7^ are new arbitrary constants. 
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Thus, the first nontrivial solution of extended Gibbons-Tsarev system (35)-(37), (44), 
(45) given by (69), (74) and (75) leads to the first nontrivial component hydrodynamic 
reduction (57), (58) of the Manakov-Santini system. In such a case, the generating 
function of conservation law densities can be found explicitly 

M M N-M N-M 

Inp = - 5 ^ e^(a- - g) - ^ 7^ ln(a"^ -q)-qY,Pn ln(&" - - ^ ^ni^^ " ^) 

m=l m=l n=l n=l 

^ M N-M 

+ i E E - Q) + 21n(a™ - q) ■ ln(6" - g) - In^a^^ - q)] (76) 

m=l n=l 

- In 1 + y + - y y f^n^m [ Li2^-^ - n^-^ — ^ ) , 

\ m=\ ^/ m=l n=l ^ ^ 

where the bi-logarithm 

OO J. 

k=\ 

A corresponding expression for the function F is given by 

^ M M N-M N-M 

F = -q' + 6j2Ua'^-Q) + Yl - + ? E " + E ^-(^^ " ^) 

m=l m=l n=l n=l 

M N-M 

+ i E E (^nem[ln\a^ - g) - 2 ln(a- - q) ■ ln(6" - g) - ln2(&" - g)] (77) 

m=l n=l 

+ 2 E E /^n^- (Li,-^ - Li,^ j . 

m=l n=l ^ ' 

9 Generalized hodograph method 

In comparison with integrable hydrodynamic chains their semi-Hamiltonian hydrody- 
namic reductions possess infinite series of conservation laws and commuting flows (cf. 
(24) and (28)). Let us remind, that hydrodynamic type system (57), (58) is associated 
with the Riemann surface A = A(a; g), whose M branch points A* = A(a; g*) determined 
by M solutions g*(a) of the algebraic equation 9A(a;g)/9g = are first M Riemann 
invariants; all other N — M Riemann invariants = A (A; IP) are nothing else but just 
mark points on this Riemann surface. 
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Since a dependence a, b) can be inverted to a, b) for any hydrodynamic re- 
duction (57), (58), generating function of conservation laws (26) 

dtP = d.x[{q{p) - a)p] 

leads to infinite series of conservation laws (see [30]) by virtue of formal expansions 

p^'^ = p^+p;i)A«+p;2)(A«)'+p;3)(A^'^)' + ..., ^ = l,2,...,M, 

(78) 

= p^„+pJ-^^A^^Vpj2)(A^'^)'+p;3)(^^'^)' + -- J = l,2,...,iV-M, 

where A*-*^ and A^'^^ are local parameters at the vicinities of Po^ respectively. How- 

ever, in such a case, dependencies qiPo) as well as q{po) are highly complicated, because 
the dependence p{q) is much more simpler than q{p) in all known examples (see the 
previous Section). Thus, instead (78), we utilize another formal expansions given by 

g» = a^ + a;i)A(^) + a;2)(A(^))2 + a^(3)(A»)3 + ..., z = l,2,...,M, 

(79) 

gO) = 6.+6;-^^A^^V6j2)(A^''^)' + 5;3)(^^'^)' + -> J = l,2,...,iV-M, 

where A^*^ and a''"'"' are local parameters at the vicinities of a* and 6^ , respectively. Since, 
the dependence p{q) possesses M singularities in points a* and N — M singularities in 
points 6^, formal expansions (78) cannot exist in a general case. Nevertheless, in some 
special cases, first N conservation laws can be found in the form 

dtPl = d,[{a' - a)pi] , dtpl = d,[{b^ - a)fP,] . (80) 

Below we restrict our consideration on a simplest sub-case (67), which is compatible 
with any reduction of X{p) (see the previous Section). In this case, infinite series of 
conservation laws can be found utilizing (79). It means that first conservation laws are 
determined by (80). 
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Without loss of generality, let us choose the equation for Riemann surface (70). In 
such a case (see (67) and (69)), three dimensional four component hydrodynamic type 
system (7) reduces to the pair of commuting component hydrodynamic type systems 
(see (9) and (57), (58)) 

a\ = (a* — a)a^ + Ma:, Qy = {{dY — aa^ — c)al. + a^Ux + Ut, i = l,2, ...,M, 

Iji = (^ifi - a)bi + Ux, Vy = {{Vf -aV -c)V, + Vux + Uu j = 1,2, N - M, 
where (e = Se^, /3 = 7 = E7^) 

^ M M N-M 

m=l m=l n=l 

Then all coefficients p*^^^ and can be found iteratively, while coefficients a*^^j(a, b) 
and 6|^.|(a, b) are determined by formal expansions of equation A = A(a; g) of a Rie- 
mann surface. A substitution of these expansions g*^'^ and q^^'^ together with other four 
expansions 

= dto. + \'^''>dti,. + {\'^'^fdt2, + {\^'^fdt3,^ + ..., 
5« = .io)+^;i)^^^^ + ^(2)(^^^'^)' + ^(3)(A«)' + -, 

= 5;o) + 5j,)A^^V.l,)(A^^'V + §j3)(A^^'V + ... 

to (24) leads to A^ infinite series of commuting flows. Then the generalized hodograph 
method admits to construct inflnitely many particular solutions in an implicit form (see 
[35]). 
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10 Conclusion 



Manakov-Santini system (7) is the first example in the theory of three dimensional hy- 
drodynamic type systems, which is naturally equipped by the two component pseudo- 
differentials g, /, determined by systems in partial derivatives of the first order 



(q-a) 



u 

+ ' 



f 



(g^ -aq-c) 



( \ 



I \ 

u 



We believe that more complicated vector pseudo-differentials 



/i = (?(/, u)/,. + F(/,u)n, 



/ \ 

u 



are associated with vector hydrodynamic chains (4). 
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